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INTRODUCTION 


This  document  is  the  Summary  Progress  Report  submitted  by  Robert  Gold 
Associates  in  accordance  with  the  CDRL  of  Contract  No.  N001 73-80-C-0480  and 
describes  the  work  accomplished  over  the  reporting  period  of  September  17,  1980 
through  July  1,  1981  on  the  Study  of  Multistate  Pseudonoise  Sequences. 

Although  mainly  binary  encoding  sequences  have  been  used  to  date  as 
the  encoding  mechanism  for  spread  spectrum  systems,  and  their  employment  has 
resulted  in  striking  improvements  in  performance  and  operational  capabilities 
over  conventional  communica tion  systems,  it  has  been  recognized  that  further 
substantial  improvements  in  spread  spectrum  communication  systems  could  be 
economically  achieved  through  an  increased  knowledge  and  application  of  multi¬ 
state  (as  distinguished  from  binary)  sequences. 

-  The  present  work  has  been  successful  in  extending  the  fundamental  con¬ 
cepts  and  results  of  the  theory  of  binary  sequences  to  the  multistate  case. 
Substantial  progress  and  new  results  have  been  achieved  in  the  following  areas: 

1.  Generation  of  Maximal  Linear  Multistate  Sequences 

2.  Implementation  of  Multistate  Sequences  with  Binary  Logic 
Elements 

3.  Special  Properties  of  Multistate  Sequences 

4.  Definition  and  Computation  of  the  Autocorrelation  Function 
of  Multistate  Sequences 

5.  Computation  of  the  Power  Spectral  Density  of  Maximal  Linear 
Multistate  Sequences  • 

In  Section  1  the  multistate  encoded  signal  and  its  analytical  repre¬ 
sentation  is  reviewed  to  ensure  that  the  properties  of  the  multistate  sequences 
which  have  the  greatest  impact  on  the  resultant  signals  were  emphasized  during 
the  current  study. 

In  Section  2  the  concept  of  a  linear  multistate  sequence  is  discussed. 
The  multistate  alphabet  (finite  fields  over  (IF (2  ))  and  the  polynomial  corres¬ 
ponding  to  a  feedback  shift  register  are  described.  A  method  for  implementing 
multistate  sequences  with  binary  logic  elements  is  presented  in  this  section. 
This  result  assures  the  practical  and  economical  application  of  multistate 
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sequences.  The  concept  of  a  maximal  multistate  sequence  is  discussed,  and  a 
method  for  determining  maximal  multistate  polynomials  is  given.  A  table  of 
irreducible  polynomials  of  degree  2  for  the  8-state  case  is  constructed  and 
presented  in  Section  2.5. 

The  autocorrelation  function  of  multistate  sequences  is  defined  in 
Section  3,  and  the  autocorrelation  function  of  maximal  linear  multistate 
sequences  is  determined.  The  result  is  essentially  independent  of  the  period 
of  the  sequence  and  depends  only  on  the  number  of  states  of  the  sequence.  It 
is  found  that  the  autocorrelation  function  of  multistate  sequences  has  regu¬ 
larly  spaced,  relatively' high  out-of-phase  autocorrelation  peaks.  In  Sec¬ 
tion  3.8  we  describe  how  these  out-of-phase  correlation  sidelobes  of  a  multi¬ 
state  maximal  linear  sequence  may  be  used  as  an  aid  to  rapid  synchronization 
of  pseudonoise  encoded  systems. 

The  power  spectral  density  of  maximal  linear  multistate  sequences  is 
determined  in  Section  3.8  where  a  formula  for  the  spectral  lines  is  given  in 
terms  of  the  discrete  Fourier  transform  of  a  sequence  of  2  -  1  values,  where 
2^  is  the  number  of  states.  This  formula  allows  the  easy  computation  of  the 
spectrum  of  multistate  encoded  signals.  The  power  spectral  density  of  a  maxi¬ 
mal  8-state  sequence  of  period  63  is  computed  and  presented  in  this  section. 
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1.  DEFINITION  OF  MOLT  1  STATE  SIGNAL  AND  ITS  CORRELATION  FUNCTION 

Although  tlie  current  study  is  concerned  with  the  properties  of  the 
digital  multistate  sequences,  which  are  the  Sasic  encoding  mechanism  for 
spread  spectrum  signals,  it  is  important  to  define  the  role  of  these  multi¬ 
state  sequences  in  determining  the  properties  of  the  resultant  signals.  This 
will  ensure  that  the  direction  of  the  current  work  can  give  priority  to  the 
study  of  those  properties  of  the  multistate  sequences  which  have  the  greatest 
impact  on  the  resultant  signals  and  their  application  to  tactical  communica¬ 
tion  systems.  Accordingly,  a  review  of  the  analytical  representation  of 
pseudonoise  signals  and  their  correlation  function  was  undertaken  and  docu¬ 
mented.  The  basic  formulas  are  summarized  in  what  follows. 


A  pseudonoisc  signal  which  has  been  encoded  by  a 
multi  state  pseudorandom  sequence  may  be  represented 
analytically  as 

iV* 

s(t)  =  2  COS  (lot  +  Y(a(k)))  f0(t-Ak) 

k  =  -  » 


Re 


i  wt 


2  e 

k=  -  oo 


1  yC3 (k) 3  fQ ( t  _  Ak) 


where  m  =  carrier  frequence 

a  =  multi st ate  sequence 

Y  =  real  valued  mapping  of  multistatc  sequence 
which  determines  how  the  multistate  sequence 
will  operate  on  the  carrier 

f  =  squarewave  function;  f()ft)  =  1;  0  <  t  <  A 

f  ()  (  t  1  =  0  ;  otherwise 

A  =  code  chip  width. 
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Example:  Let  a  be  the  4-state  sequence  with  values 

0,  1,  6,  B2. 

Let  y  be  the  mapping  such  that 
Y  : 


B 


0 

1 

B 

2 


7T 

4 


-  t r 

X 

3tt 

4 


IT 


The  resultant  phase  modulated  signal  is  illustrated 
in  Fi gure  1 . 


Figure  i.  Four-phase  modulated  carrier 

The  correlation  function  of  wideband  PN  signals  is 
a  crucial  parameter  in  measuring  their  performance  in 
spread  spectrum  systems.  The  cros s - cor re  1  at i on  function 


of  two  pseudonoisc  signals  s,(t)  and  s,(t)  which  have  been 

n  d 

encoded  by  two  multistate  sequences  a  and  b,  respectively, 
has  been  shown  to  be  given  by 

oo  j  .  P  "  1  * _ 

R(sa,sb)(T)  =  ^  4Re|cla,Ti^^n(a(k))n(b(k-il))jR(f0,f0)CT-Ak) 

where 

n(a)(k)  =  and  n(b)  = 

in  view  of  this  result,  the  correlation  function  of 
two  multistate  sequences  has  been  defined  as  follows  [1, 
p.  97]. 

Definition:  Let  a  and  b  be  two  multistate  sequences 

of  period  p.  Let  n  be  any  complex  valued  mapping  of  the 
sequence  alphabet.  The  correlation  function  of  the  multi - 
state  sequence  with  respect  to  the  mapping  n  is  defined 
to  be 

t  P- 1  - 

0 n  ( a  , b  )  ( r )  =  i.  y  n  (a(k) )  n  (b  (k-  t)  }  • 
n  p  k  =  0 

Hi  us  , 

CO 

R(sa,sb]  =  y  i-Re{olwT0  (a,h)(e)}R(f  ,f  )(r-Al) 

I  =  -  “  * 


1  -  1  Detection  of  Signal  Encoded  with  Multistate  Signal 

The  received  signal  is  given  by 


s 

a 


(t) 


E  cos/w^t  +  Y(a(b)) 
k=-^>  V 


f0(t  -  Ak) 


i  w  t  ,r' 

c  E 


iV(a (k) ) 


1 


e 


e 


•  f^it  -  Ak) 


Let  the  locally  generated  signal  be 


Sb(t-T)  = 


i  cos(( 


(^c  —  oj  i )  t  +  y  (b  (k) ) )  fn(t-x 


>)  V 


where  w  ^  is  an  IF  frequency 


y  Re 


i  (w  -  w, )  t 
e  L  e 

k=-°° 


iY(b(k)) 


fo(t 


We  first  multiply  the  received  signal  s  (t)  by 
generated  signal  s^(t-  -  t) 


s  (O 
a 


> - 


s  (t)  s  (t 
a  b 


sb(t-T) 


Le  t 


iY  (a(k) ) 

Ua(t  )  =  I)  e  fn ( t  -  Ak) 

k=-w 


0 


Ufe(t)  =  12 

k=-°° 


iy (b(k) ) 


f Q ( t  -  Ak) 


Then  we  have 


sa  ( t  >  =  2  R° 


iw  t 

e  C  U  (t) 
a 


sb(t) 


Re 


i(w  -  at  )  t 
e  C  Ub(t) 


-  Ak) 

t  -  Ak) 

the  locally 

-  t) 


s  (t)s, (t-x) 


.  iw  t  i(ui  -  u>.  )  t 

7  Re  e  C  U  (t)  Re  e  C  1 


U,  (t  -  t) 


Using  the  relationship 


Re  (  z  1  )  Re  (  z2  ) 


Z2> 


+  Y  Re(z!  '  V 


where  z^  and  z2  are  complex  numbers  and  bar  indicates  complex 
conjugation,  we  have 


Sa{t)sb(t-T> 


8  RC 


i  ( 2u)  -  to  i  )  t 

c  i 


U  (  t)  U  ( t-  T) 
a  b 


+  -Re 


U  (t)U,  (t-T) 
a  o 


We  now  bandpass  filter  about  the  IF  frequency  u) 


to  obtain 


y  (t) 


j_ 

8 


Re 


i  u 
e 


t 

U  (  t)U.  (  t  —  T  ) 
a  b 


and  bandpass  liLtor  y  (  t )  witii  a  filter  whose  lowpass  equivalent 
computes 


AP 

/  Ua(t)Ub(t-T)dt 

0 


Ibis  integral  has  been  shown  to  be  equal  to 


T:  e(a,b)(i)  0  (  f  0  ,  f0)  (t-AI) 

J  —  —  m 

For  x  =  ilk,  the  above  sum  becomes  0(a,b)(k)  and  the  output  of 
the  last  bandpass  filter  is 


7  (  t ) 


8 


R  e 


t 

9 (a , b ) (k) 


T o  envelope  detect  the  waveform  z ( t )  we  square  to  obtain 


,  1 
128  , 

Ro  J 

i  2u  t  2 

e  9  (  a  ,  b  )  (  k  ) 

+  Re 

0  (  a  ,  b  )  (  k  ) 

2~] 

< 

i  1 

J 

and  lowpass  f  i  1  t  e  r  to  obtain  -1  a->  ^J-LH  )  . . 


2 . 2  Multistat c  Feedback  Shift  Registers 

Multistate  sequences  generated  hv  feedback  shift 
registers  may  be  represented  by  polynomials  with  coeffi¬ 
cients  in  the  alphabet  field  of  the  sequence.  Examples 
are  given  below  for  the  binary,  4-state,  and  8-state  cases 

B  i  n  a  r  y 


8  -  S  t  a  t  e 


AAAAA 


a 


1 

f(x) 


6 

2 

1  +  x  +  gx" 


=  A  +  / b\  x  +  Ax2  +  Ax3  + 

The  polynomial  corresponding  to  a  multistate  shift 
register  determines  the  properties  of  the  multistate 
sequence  generated  by  that  shift  register.  Thus,  to  find 
the  period  of  any  polynomial  f  over  G!;(2n),  multiply  the 
polynomial  by  its  conjugate  polynomials  to  obtain  a  poly¬ 
nomial  f'  over  GF(2).  The  period  of  f  is  equal  to  the 
period  of  f'. 

bxamp  1  e _ (  4  -  s  t  a  t  c )  :  Gbf22)  f(x)  =  x“  +  ax  +  1 

2  7 

The  conjugate  polynomial  is  x  +  a  x  +  1 


7  2  2  4  12 

(x“  +  ax  +l)(x  +  a  x  +  1 )  =  x  +x  +  x  + 


4  .1 

I  lie  period  o  I  x  +  x  + 
■>  ? 
of  x "  +  ax  +  1  (and  x  “  + 


1 

x"  +  x  +  1  is  5  and  hence  the  period 

? 

a"x  +  lj  is  5 ,  e . g . , 


T 
X  “ 


+ 


The 


ax  +  1 


bxamp  l  e 


con  i  tig  a  t  e 


2  3  S  0 

+  ax  +  ax*”  +  x'  +  x'  +  ax  +  =  laalO  |  laalO  | 

(  8  -  s  t  a  t  c )  :  Gb  ( 23)  f  ( x )  =  x2  +  a5x  +  1 

2  3  2  6 

po  1  vnom  i  a  1  s  are  (x“  +  a  x  +  1 )  and  (x  +  a  x  +  1 )  . 
+  1  )  (  x  “  *  a'1  x  +  1  j  ( X  “  +  a(’  x  +  1  )  =  xh  +  X 3  +  1  . 
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The  period  of  x  +  x'  +  1  is  9  and  hence  the  period  of 
?  5 

f(x)  =  x“  +  u  x  +  1  (and  its  conjugates)  is  also  9,  e.g., 
1 

2  7  5  7  , 
x  +  a  x  +  1 

=  |  1  a5  a2  «  a  a2  a5  1  0  1 1  a5 
U - period  9 - fcj 


5  f  2  7  3  4  ?  s  56 

=  1  +  a  x  +  a  x“  +  ax  +ax  +  a  x  +  a  x 


+  x7  +  x9  +  a5x10 


2 . 3  Polynomials  Corrospund i to  M  n  I  L  i  s  t  a  t  e  Shift  Registers 

For  the  case  of  linear  binary  sequences,  extensive 
tables  of  irreducible  polynomials  have  been  compiled  [2] 
which  determine  the  feedback  taps  for  the  shift  registers 
which  generate  the  sequencer..  1'he  properties  of  the 
sequences  generated  by  t he  shift  registers  may  be  deter¬ 
mined  from  the  algebraic  properties  of  the  corresponding 
polynomial.  A  method  of  determining  the  polynomials  for 
multistatc  sequences  having  prescribed  properties  has 
been  developed  and  is  described  in  this  section.  The 
method  assumes  that  tables  of  polynomials  for  the  binary 
case  arc  available. 

R c s u 1 t  :  To  find  a  polynomial  f(x)  of  degree  k 
corresponding  to  a  k-stage  shift  register  which  will 
generate  a  2m- state  sequence  of  period  p,  we  proceed 
as  follows: 

(1)  Assume  a  is  a  root  of  a  maximal  binary  polynomial 
g(x)  of  degree  k*m. 

(2)  Choose  an  integer  j  such  that 

7m  •  k  j 

^  ~~  rr  ->m  •  k  77 

ged ( j  ,2  -1) 

k  -  1  i  *m 

(3)  f(x)  =  TT  (x  +  aJ  ). 

i  =  0 
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Hxamp  1c  (4-state):  in  =  2 


Wo  construct  a  polynomial  of  degree  2  corresponding 
to  a  2-stage  shift  register  which  will  generate  a  2 “  =  4 
state  sequence  of  period  p=  5. 


(1)  Let  u  he  a  root  of  the  maximal  binary  polynomial 

4  4 

of  degree  m  k  =  4  given  by  g ( x)  =  x  +  x  +  1 ,  l .  e  .  ,  a  =  a  +  1  . 


(2)  Choose  j  =  3  such  that 


,in*k  ^  -7 4  _  j 

gcd(  j  ,2m'k  -  1)  ged  ( 3 , 24 -  l ) 


(3)  f  (x  )  =  (x*a^)(x  +  a^)  =  x‘'+(a3  +  a1“)x+l 


-7  7  1  |Q 

=  x“  +  («“  +  <*  +  l)x+l  =  x"  +  a  x  +  1 


7  T  $ 

1  (x)  =  x“  +  6”x  +  1 ,  where  b  =  a‘  . 


The  shift  register  corresponding  to  this  polynomial 


The  sequences  generated  by  this  shift  register  have 
period  3 ,  c . g . , 


I  1  1  6  0  6 1  1  1 
|  period  5  | 

2 . 4  implementation  o  f  Multistate  Sequences 

[n  this  section,  we  explain  how  multistate  sequences 
may  be  implemented  with  binary  sequences  generators.  We 
shall  use  an  8-state  sequence  generator  to  illustrate  the 
general  procedure.  We  consider  the  maximal  polynomial 
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over  GF(8)  given  by  f(x)  =  ax  +  x  +  l.  The  shift  register 
corresponding  to  this  polynomial  is 


1  1 


2  2 

and  the  sequence  l/(a  x  +x+  1)  which  is  generated  by 
this  shift  register  with  the  indicated  initial  conditions 
is  given  by 


ax”  +  x  +  1 

The  complete  sequence  and  its  binary  representation  is 
shown  in  Figure  2. 

In  order  to  find  the  representation  for  the  8-state 
sequence  as  three  binary  sequences,  we  proceed  with  the 
following  steps. 

Step  1:  Rationalize  the  polynomial  representation 
of  the  8-state  sequence,  e.g., 

1  _  1  (a^x“  +  x+l)(oi^x“  +  x+l) 

2  2  *  7  2  4  1 

ax4  +  x  +  1  ax  +  x  +  1  (a“x  +  x  +  1 ) (a  x“  +  x  +  1) 

6  4  S3,  2  ,  , 

ax  +ax  +ax  +1 

6  4  3  . 

x  +  x  +  x  +  x  +  1 


The  denominator  of  the  rationalized  fraction  provides  the 
feedback  taps  for  the  required  binary  shift  register. 


r 


Figure  2.  Mnxim.il  8-State  Sequence  and  Its  Binary  Representation 


Stop  2 :  Express  the  coefficients  of  the  numerator 
from  GF(8)  as  binary  triples,  e.g., 

a6  =  [0  1  1  ]  ct  =  [010] 

a5  =  [ 1  1  0 ]  1  =  [ 1  0  0  ] 

a6x4  +  '.x5xJ  +  ax2  +  1  =  [Pll]x4  +  [110]xJ  +  [0 1 0  ]x2  +  [100] 

6  3  3  ”  ,  6  3  3  , 

x+x+x+x+1  x+x+x+x+1 


Step  3 :  Write  the  fractional  representation  of  the 
sequence  given  in  Step  2  ns  the  sum  of  three  binary  frac 
tions ,  e.g. , 

[011  ]x4  +  [1  10]x^  +  [()  10]x2  +  [100] 

0  4  3  . 

x  +  x  +  x  +  x  +  1 


x~*  +  1 _ 

xb  +  X4  +  X3  +  X  +  1 


=  1111110  1 


4  3 

X  +  X  +  X 


0  4  3  , 

X  +  X  +  X  +  X  +  1 


=  00101010 


0  4  3  . 

x  +  x  +  x  +  x+  l 


=  0  0  0  0  1  1  1  0 


.  5  ,  4  3  6  , 

11a  laaa  1 


The  implementation  of  this  8-state  sequence  generator 
using  a  simple  binary  shift  register  is  illustrated  in 


Figure  3.  Implementation  of  an  8-statc  generator  with 
a  single  binary  feedback  shift  register 

2.5  Construction  of  Tables  of  Irreducible  Polynomials  of 
Degree  2  over  G  F  (  2  '■* ) _ 

In  order  to  illustrate  the  computation  of  useful  tables 

of  multistate  irreducible  polynomials,  and  in  order  to  have  at 

hand  a  convenient  list  of  such  polynomials  to  aid  in  tiie  study 

of  multistate  sequences,  a  table  of  irreducible  polynomials 

3 

of  degree  2  over  the  finite  field  GF(2  )  was  constructed. 

These  polynomials  will  generate  eight  (8)  state  pseudonoise 
sequences  whose  periods  are  factors  of  63. 

The  multistate  polynomials  in  each  group  of  three  (3)  of 
Table  1  are  factors  of  the  binary  polynomial  in  the  product 
column.  T  h u  s ,  for  example 

(x2  +  p.3x  +  B)(x2  +  p/’x  +  P2)(x"  +  P3x  +  B2)  =  x6  +  x  +  1 
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Table  1.  Irreducible  Polynomials  of  Degree  2  over  GF(2  ) 


Period 


Product 


2 

3 

1 

X 

+ 

6  x  +  P 

2 

6 

2 

X 

+ 

6  x  +  P 

2 

5 

4 

X 

+ 

6  x  +  6 

2 

5 

3 

X 

+ 

e  x  +  e 

2 

3 

6 

X 

+ 

6  x  +  6 

12 

2 

X 

+ 

P6x  +  6 

2 

2 

5 

X 

+ 

ex  +  e 

10 

2 

X 

+ 

e 4  x  +  e 

20 

2 

x 

+ 

Px  +  e6 

7 

14 

28 


5 

x  +  e  x  +  l 


X4  +  ?  X  +  1 


x 2  +  £  ^  x  +  1 


63 


21 


63 


x  +  x  +  1 


6  ^  4  2  ,  .  . 

X  +  X  +  X  +  X+1 


6  ,  5  ,  2  ,  ,  , 

X  +  X  +  X  +  X+1 


X6  +  X3  +  1 


11 

22 

44 


x2  +  x  +  g4 


x  +  x  +  6 

2  2 
x  +  x  +  f 


63 


6  5  3,  2  . 

x  +  x  +  x  +  x  +1 


1  3 

26 

52 


x2  +  e6x  +  e6 
x2  +  e5x  +  e3 
x2  +  e3x  +  e3 


63 


6  .  5  3  2  . 

X  +  X  +x  +  x  +1 


15 

30 

60 


x  +  6  x  +  B 

2  2 

x  +  ex  +  g 

2  2  4 

x^  +  e  x  +  e 


21 


6  A  5  4  2  . 

X  +  x  +  x  +  x  +1 


21 


X  +  X  +  1 


X  +  X  +  1 
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Table  1.  Irreducible  Polynomials  of  Degree 


I 


2 

4  ' 

23 

X 

+ 

6  x  +  6 

46 

2 

X 

+ 

4 

6x  +  6 

2 

.2 

29 

X 

+ 

6  x  + 

2 

„  .  3 

31 

X 

+ 

6x  +  6 

2 

62 

X 

+ 

6“x  +  6 

61 

2 

X 

+ 

4 

e  x  +  e 

63 


63 
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The  roots  of  the  polynomial  x  +  x  +  1  are 

,  2  A  8  16  32, 

<  a  a  a  a  a  a  J 


and  a 

8 

of 

2 

3 

are 

root  s 

x 

+ 

£ 

+ 

B 

and 

16 

of 

2 

6 

„  2 

9 

a 

are 

roots 

X 

+ 

B 

+ 

C 

B  =  a 

and 

32 

of 

2 

5 

U 

3  „  2 

a 

are 

roots 

X 

+ 

£ 

+ 

a 

£  +  £  +  1  =  0 

2.6  Some  Properties  of  Multistate  Sequences 

2.6.1  Multiplication  of  a  Multistate  Sequence  by  a  Character 

The  following  result  describes  what  happens  when  a  multi-state  sequence 
is  multiplied  by  a  character.  The  resulting  sequence  is  a  phase  shift  of 
the  original  sequence.  The  precise  phase  shift  obtained  is  described  by 
the  following  result. 


Result:  Suppose  we  have  a  2  -state  maximal  sequence  h  corresponding  to 

the  nth  degree  polynomial  f(x)  e  GF(2k)[x]  with  root  aV  in  GF(2kn). 

2k  -  1 

(The  values  taken  by  this  sequence  are  {0,  1,  6  ...  £  ) 

2k‘n  -  l  M  i  -1 

where  £  =  a  — - -  =  a  .)  Then  8J  h  ( i )  =  h(i  -  -jv  M)  . 

2-1 


Proof:  By  previous  result  there  exists  a  linear  mapping 

k  k  ~i 

if  of  GF(2  -n)  into  GF(2  )  such  that  h(i)  =  f((a%)  ). 


_ 

1 

-  jv 

Mj  = 

-(i 

-  jv 

aV) 

-i  *1 

aV) 

,  Ml 

(a  )! 

-i 

n 

aV) 

£J  = 

=  £J  h  ( i ) 


Example:  (see  following  page) 


19 


2.6.2  Binary  Representation  of  Maximal  Multistatc  Sequences 

3 

Result:  Suppose  we  have  a  2  -state  maximal  sequence  h  corresponding  to 

the  degree  polynomial  f(x)  c  GF(2^)[x]  with  root  aV  in  GF(23n). 
g  (x) 

Let  be  a  rationalized  representation  of  h.  Then  the  binary 

representation  of  the  sequence  h  =  is  given  by 


h  =  M-  pn  ,  *  pih 

'  7  7  —  ( v  M)  '  7 


(v  M) 


where:  M  = 


23n  -  1 
23  -  1 


Proof :  h  = 


'Pj(g(x))  P2(g(x))  P3(g(x))] 


[Pjigfx))  Pj  (£_1(g(x)  ))  PjC6(g(x))) 


0}  pi((e_1>(ffS))- 


p.  ,f)  frS 


*b)-(F'h),-m-  (p>h) 

-(v  M) 


Example:  Suppose  we  consider  the  2  -state  maximal  sequence  h  correspond¬ 

ing  to  the  2n3  degree  polynomial  x"  +  a^x  4  a^.  The  sequence  gener¬ 
ated  by  the  polynomial  and  the  binary  representation  of  the  sequence 
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have  been  shown  in  Section  2.4  to  be  the  sequences  of  Figure  2.  We  note 
from  Table  1  that  a*3  is  a  root  of  the  polynomial  x2  +  a^x  *-  .  Thus 

v  =  13;  v  =34  and  v  3M  =  54.  Thus  the  binary  sequence  representation 
of  the  8-state  sequence  is 


This  may  be  verified  by  observing  the  binary  sequence  representation  of 
Figure  2. 

2.6.3  Sequences  Generated  by  Conjugate  Polynomials 

R 

The  k-1  conjugate  polynomial  to  a  polynomial  f(x)  e  GF(2  )  are  obtained 

2 

by  replacing  the  coefficient  f(j)  of  x^  by  (f  ( j  ))  *"  ,  (f(j))2  ,  .  .  . 

2k_1 

f(j)  *■  ,  respectively.  Thus,  for  example,  the  conjugate  polynomials  to 

x2  +  B^x  +  B^  c  GF(23)  [x]  are  found  to  be 


l2x2+  (B6)  x+  (fV°  -  x2  +  B5x  +  B5 


l4  X2  +  (fS  x  +  (f6)  =  X2  +  e3  +  B3 


The  2  state  sequences  generated  by  these  conjugate  polynomials  can  be 

obtained  from  the  sequences  h  generated  by  the  polynomial  f(x)  by  replacing 

2  k-1 

each  term  h(i)  of  the  sequence  by  (h(i))~,  (h(i))~  ,  .  .  .  (h(i))2  , 

respectively.  Thus  in  the  above  example  if  h(0),  h(l),  ...  is  the 

sequence  generated  by  f(x)  =  x“  +  u^x  +  ,  then  (h(0))",  (h  ( 1 ))  .  .  .  and 

(h(0))4,  (lid))4,  .  .  .  will  be  the  sequences  generated  bv  x~  +  a3x  +  a3  and 
o  3  3 

x  +  a  x  +  n  ,  respectively. 
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Since  the  j1*1  conjugate  polynomial  to  the  polynomial  f(>:)  f  C.F(2  ) 

,  j  t 

having  root  aX  has  root  v,  the  binary  representation  of  the  sequence  cor¬ 
responding  to  this  conjugate  polynomial  is 


^66  13 

Example:  The  polynomial  x  +  £  x  +  £  having  root  a  will  generate  the 

5  A  3 

sequence  h  =  1  1  £  1  £  £  £  ....  This  sequence  can  be  repre¬ 

sented  by  phase  shifts  of  a  binary  sequence  a  corresponding  to  the 


polynomial  x*3  +  x-"  +  x  3  +  x2  +  1  =  (x2  +  £  x+£  )  (x  +  £3x+B  )  (x  +B  x+B3). 


In  fact  in  Figure  2  we  have 


h 

=  i 

1 

„5 

D 

i  eA 

e3 

a 

=  i 

1 

1 

i  i 

i 

a9 

=  0 

0 

1 

0  1 

0 

a  „ 

=  0 

0 

0 

0  1 

i 

-y 

The  conjugate 

polynomial 

x2+£5x 

f  £6 

sequence 

1 

h  = 

1  1 

B6  1  F 

£6  £5  . 

•  •  9 

sequence 

h  by 

replacing 

each  term  h(i 

0  ... 

1  ... 

1  ... 

2  6 

having  root  a  will  generate  the 
which  is  obtained  from  the 
by  (Vi(i)j".  The  sequence  h 
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can  be  represented  by  phase  shifts  of  the  sequence  a  in  accordance  with 
the  above  result.  Thus  we  have 


? 

h 

1 

1  f3 

1 

e  e6 

E5 

.  .  . 

a  1 8 

t 

a  = 

1 

1  1 

1 

0  0 

1 

.  .  . 

a32x9 

a3fc  = 

0 

0  0 

0 

l  i 

1 

.  .  . 

a- 32x9 

a-36  = 

0 

0  1 

0 

0  1 

0 

.  .  . 

The  conjugate 

polynomial 

2 

X 

3  3 

+  ex  +  e 

having  root 

52 

a  ’ 

will  generate 

the  sequence  h 

t  1 

=  116 

6  1 

2  5  3 

e  e  b  . 

•  • 

.  which  is 

obtained  from  the 

sequence  h  by 

replacing  each 

term  h(i) 

by 

(h(i))\ 

The 

1  F 

sequence  h 

can  be  represented  by  phase 

shifts  of 

the 

sequence 

a  in 

accordance 

with  the  above 

resul t . 

Thus 

we  have 

t  t 

h 

1 

i  e6 

1 

b2  e5 

e3 

.  .  . 

a27  = 

1  f 

a  = 

1 

1  0 

1 

0  1 

i 

.  .  . 

a16x9  = 

t  i 

al  8 

0 

0  1 

0 

0  1 

0 

.  .  . 

a-16x9  = 

«  f 

a-l  8  = 

0 

0  1 

0 

1  0 

i 

.  .  . 

v 


♦ 
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3. 


CORRELATION  FUNCTION  OF  MULTI  STATE  SEQUENCES 


V 

i 


I 


3 . 1  Definition  of  9 

The  correlation  function  with  respect  to  a  complex  valued  mapping  of  the 
sequence  alphabet  of  two  (2)  multi-state  sequences  has  been  defined  (in  Report 
No.  1)  as  follows 


6  (a  ,b)  (t ) 
n 


t’-i 

=  ^  ^  ^  n(a(k))  r(b(k-r)) 


where  a,b  are  multi-state  sequences  of  period  P 

n  is  complex  valued  mapping  of  the  sequence  alphabet 

In  what  follows,  we  compute  the  auto  correlation  function  with  respect  to 
a  specific  complex  valued  function  n  of  maximal  multi-state  PN  sequences. 

3 . 2  The  Mapping  of  n 

The  specific  complex  valued  mapping  of  CF(2  )  which  we  shall  use  has 
values  on  the  complex  unit  circle.  The  mapping  n  ultimately  determines  the 
correspondence  between  the  characters  of  the  multi-state  alphabet  and  the  phase 
of  the  carrier  which  will  be  modulated  by  the  multi-state  sequence. 

Definition:  n  is  a  complex  valued  mapping  of 

CF(2k)  =  jo,  1,  B,  E"  .  .  .  62  “  >  given  by 


n  (0)  1 


n(6s) 


We  note  that 


i  (s  +  1  )tt 
.  k-1  ' 


0  <  s  <  2  -2 


n(BS) 


i  (s  +  2)m 
~k- 1 


2k  -  1  1  s  i  2k+1  -  3 


■?  fc  ' 


6 (a ,a)  (t) 


-1 


otherwise 


2n‘k-l 


6  (.a, a)  (t)  =  n(a(t)^  n(a(t-t)j 

t=0 


definition  of  £ 


6  (a , a) (jM) 


2n ' k- 1 


^  ^  n  (a ( t ))  n(a(t-jM)J 


t=0 


2n ' k- 1 


^  ^  n  (a  ( t )) 


n  a  (t  -  ( jv) (v 


■'><)) 


t=0 

2° ' k- 1 


^  ^  n(a(t))  a(t)^ 


t=0 


M  = 


2k  n-I 
2k-l 


previous  result  /'  (3.0) 


For  each  non-zero  term  a(t)  e  GF(2  )  of  the  multi-state  sequence  a 
there  is  an  integer  f  ^a  ( t ))  where  0  2  E  (a  ( t  jj  -  2k-2.  Such  that 


(a(tl) 


a(t)  =  6  V  >  .  Thus  we  have 


2n" k-l 


y^n(sUa(t)))r,(ejv  +  t(a(t)))  +  2k 


(n-1) 


-  1 


t=0 

a(t)/0 


2n-k-l 

£ 

a(t)*0 


i(t(a(t))  +  l)  it 

7k -1 


i  (f  ( a  ( t ) )  +  v  j  +  7i 


,k-l 


jv+  JL( a  ( t )  )  <  2  -2 
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2nk-l 

e 

a(t)^0 


i 

J(a(t))  +  l)* 

i  (i  (a  ( t ) )  +  vj  +  2y 

71 

e 

2k- 1 

2k-l 

e 

+  2k(n-l)  _  j 


jv+  £.(a(t) )  >  2  -2 


2n ’ k- 1 

■  E 

a(t)/0 
i (a(t) )  1  2k— 2— j 


i  v  j  it 

2  "  - 1 

V' 

i  (vj  +  1  )ir 

2k_1 

+  X 

2k_1 

a 

o 

II 

<3 

e 

u 

+  2k(n-1}  -  1 


£ (a  ( t ) )  >  2  -2-vj 


k  k 

Since  there  are  (2  -l)-vj  values  for  which  J(a(t))  i  2  -2-vj,  and 

there  are  vj  values  for  which  £(a(t))>  2  -2-vj,  and  since  each  value 

.k(n-l)  . 

occurs  in  the  sequence  2  times,  we  have 


0 (a , a) (jM) 


(2k-‘-vj) 


k(n_I}_i^p  k(n.1}  .  JL(vj±i)l  k(n_i) 

2  i  /  . ,  2 


+  (vj)  2  e 


+  2 


3 . 4  Autocorrelation  Function  for  8-State  Maximal  Sequences 

For  8-state  maximal  sequences,  the  above  correlation  formula  becomes 


23n-l, 


23-l 

© |.i  J=  (7-  vj)  2  e 


3(n-J)  3  ( n  —  1 )  .  iMtllZ.  3(n-l) 

A  M 


+  (vj)  2  e 


+  2 


1 


For  a  two-stage  maximal  register  which  generates  an  eight-state  sequence 
of  period  63  (n=2),  the  correlation  formula  becomes 


i  v  j  t  i  (vj  +  1  )tt 

4  4 

6(9j)  =  (7-vj)8e  +  vj  8  e  +7 


Example:  We  use  the  above  formula  to  compute  the  auto  correlation  function 

of  the  maximal  8-state  pseudo  noise  sequence  corresponding  to  the 
polynomial  x3  +  a^x  +  a3.  We  note  from  Table  2-1  that  a33  is  a  root 
of  this  polynomial.  We  thus  have 


28 


i  2  i  ti 

(21+  l)w 

e(9j)  = 

(7  -  2 j )  (8)  e  4  + 

(2j)(8)  e 

4  +  7 

1271 

1  3t 

e  (9)  = 

(5) (8) 

e 

4  +  (16)  e 

4  + 

7  = 

(7  -  8/2)  +  i (-40  -  8/2) 

i4r 

i5*n 

6(18)  = 

(3) (8) 

e 

4  +  (32)  e 

4  + 

7  == 

(-17  -  16/2)  +  i  1  bfi 

i6r 

i  7  TT 

e (27)  = 

( 1) (8) 

e 

4  +  (6) (8)  e 

4  + 

7  = 

(7  +  24>/2)  +  i  (8  +  24. 2) 

ir 

i  2rr 

6(36)  = 

(6) (8) 

e 

4  +  ( 1 ) (8)  e 

4  + 

7  = 

(7  +  24/2)  +  i  (-8  -  24/2) 

1  3tt 

i4r 

6(45)  = 

(4) (8) 

e 

4  +  (3) (8)  e 

4  + 

7  = 

(-17  -  lb.'T)  +  i(-16./2) 

i5r 

6  r 

6(54)  = 

(2) (8) 

e 

4  +  (5) (8)  e 

4  + 

7  = 

(7  -  8.2)  +  i  (40  +  8.1) 

Example:  We  use  the  above  formula  to  compute  the  auto  correlation  function 

of  the  maximal  8-state  pseudo  noise  sequence  corresponding  to  the 

2  4  2  23 

polynomial  x  +  6  x  +  £  .  We  note  from  Table  1  that  a  is  a  root 

of  this  polynomial.  We  thus  have 

_  i2i~  _  i  (2.i  +  1  )t 


e  (9  j ) 

e  (9) 

6(18) 

6(27) 


(7  -  2 j ) (8)  e  +  (2j ) (8)  e 

i  2  it  i  3tt 

(5) (8)  e  4  +  16  e  4  +  7  = 

i4r  i 5r 

(3) (8)  e  4  +  32  e  4  +  7  = 

ifer  i  7  tt 

( 1 ) (8)  e  4  +  (6)(8)  e  4  +  7  = 


+  7 

(7  -  8^)  -t  i  ( -40  -  8/2)  =  (-4  •  31  -  i  (5 1  •  31  j) 
(-17  - 1 6 »/2‘)  +  1(16/2  =  (-39-63+i  22  •  63) 


(7  +  2 4 *'2 )  +  i(8  +  24/2)  *=  (40  •  94  +  i  41-94) 


|  0(a,a) | 


Figure  U.  Absolute  Value  of  Autocorrelation  Function  of  Maximal 
Linear  Sequence  of  Period  63  over  GF(23) 


Sequence  Period 


Autocorrelation  Function  of  Maximal  Linear  Sequence  over  GF(2  )  (Eight  State) 


XTT  12^ 

6(36)  =  (6)  (8)  e  4  +  (1)(8)  e  4  +  7 

i  3tt  i4r 

6(45)  =  (A) (8)  e  4  +  (3) (8)  e  4  +  7 

_  i5tt  i6r 

6(54)  =  (2) (8)  e  4  +  (5) (8)  e  4  +  7 

A  plot  of  the  absolute  value  |6(a,a)|  is  presented  in  Figure  4,  and  of 
9(a,a)  in  Figure  5. 

3.5  Computation  of  the  Autocorrelation  Function  of  16-State  Maximal  Linear 
Sequences _ 

In  this  section  we  give  examples  of  the  computation  of  the  autocorrela¬ 
tion  function  of  a  16-state  maximal  linear  sequence.  The  unnormalized  auto- 

k  k  •  n 

correlation  function  of  a  2  -state  maximal  linear  sequence  of  period  2  -  1 

has  been  shown  to  be  given  by  the  formula: 

0(a,a)(jM)  =  2k^n'^  Qk(rj)  -  1 

9(a,a)(  t  )  =  -  1  r  £  jM 


=  (7+  24/2)  +  i(-8  -  24/2)  =  (40  •  94  -  i  41  •  94) 

=  (-17  -  16/2)  +  i(-16/2)  =  (-39  •  63  -  i  22.  •  63) 


=  (7  -  8/2)  +  i  (40  +  8/2)  =  -4  •  31  +  i  51  •  31 


where  M 


2k ' n  -  1 
2k  -  1 


Sequence  Period 


(No.  of  Sequence  States)  -  1 


No.  of 

Correlation  Peaks 


and  Q.(j; 


(2k  - 1  -  j) 


i  J 

ok  -  1 


Ki +  U17 

ok  -  1 


+  j  e 


+  1 


Q4(j)  = 


-  1  J  * 

,  8 

(15  -  j  -  l)  e  +  j  e 


jliijjz 

8 


+  1 


32 


J 


Q4(j) 


1 

14.6 

- 

i 

6.1 

2 

11.0 

- 

i 

11.0 

3 

5.6 

- 

i 

14.1 

4 

-  0.5 

- 

i 

14.7 

5 

-  6.4 

- 

i 

12.8 

6 

-10.9 

- 

i 

8.7 

7 

-13.4 

- 

i 

3.1 

8 

-13.4 

+ 

i 

3.1 

9 

-10.9 

+ 

i 

8.7 

10 

-  6.4 

+ 

i 

12.8 

11 

-  0.5 

+ 

i 

14.7 

12 

5.6 

+ 

i 

14.1 

13 

10.9 

+ 

i 

11.0 

14 

14.6 

+ 

i 

6.1 

The  above  values  of  determine  the  auto  correlation  function  of  16-state 
maximal  linear  sequences  of  any  period.  Thus,  for  example,  for  n  =  2, 
period  255  we  have: 


j 

9  ( 1 7  J ) 

16  q4(j)  -  l 

|6(17j)| 

1 9( 17 j ) 1/255 

1 

232  .6 

- 

i 

97.6 

252.2 

.99 

2 

175.0 

- 

i 

176.0 

248.2 

.97 

3 

88.6 

- 

i 

225  .6 

242.4 

.95 

4 

-  9.0 

- 

i 

235.2 

235.4 

.92 

5 

-103.4 

- 

i 

204.8 

229  4 

.90 

6 

-175.4 

- 

i 

139.2 

224.0 

.88 

7 

-215.4 

- 

i 

49.6 

221.0 

.87 

8 

-215.4 

+ 

i 

49.6 

221.0 

.87 

9 

-175.4 

+ 

i 

139.2 

224.0 

.88 

10 

-103.4 

+ 

i 

204.8 

229.4 

.90 

11 

-  9.0 

+ 

i 

235.2 

235.4 

.92 

12 

88.6 

+ 

i 

225.6 

242.4 

.95 

13 

175.0 

+ 

i 

176.0 

248.2 

.97 

14 

232.6 

+ 

i 

97.6 

252.2 

.99 

The  absolute  values  of  9  are  plotted  in  Figure  6. 
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The  expression  for  the  normalized  auto  correlation  function  is: 


9(a,a)( jM)  = 


2  KM  - 

kn 

2  -  1 


This  expression  may  be  approximated  by: 


9(a,a)(jM)  - 


QK(rj) 


We  note  that  this  expression  is  independent  of  n. 


3.6  Asymtotic  Expression  for  9^ 

The  absolute  value  of  the  normalized  auto  correlation  peaks  of  multi¬ 
state  maximal  linear  sequences  approach  1  as  the  number  of  states  increases. 
In  fact  we  have: 


Proof: 


9(jM) 

9(0) 


9(jM) 

9(0) 


->  1  for  large  k 


qk(j)' 


1  3  T  _  M,  j  +  1  lz 

k  \  ?k-l  *  ?k-l 

2  -  1  -  J  e  +  j  e  +  1 
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1  j  TT  1'  j  TT 

i  - 1 

/  k  \  J<-1  jc-l 

1 

7T- 

1 

V2  -  V  e  2  -  j  e  2 

1  -  e  2 

+  1 

3.7  Use  of  Correlation  Sidelobes  for  Signal  Acquisition 

In  this  section  we  describe  how  the  out-of-phase  correlation  sidelobes 
of  a  multi-state  maximal  linear  sequence  may  be  used  as  a  synchronization  aid. 
We  show  that  the  time  required  for  the  synchronization  of  multi -state  sequences 
is  significantly  less  than  that  required  for  the  synchronization  of  binary 
maximal  linear  sequences. 

The  conventional  synchronization  techniques  for  binary  sequences  compares 
the  received  code  with  successive  phase  shifts  of  a  locally  generated  replica 
of  the  received  code  until  phase  synchronism  of  the  two  (2)  codes  is  achieved. 
If  the  codes  have  period  N  and  the  integration  is  carried  out  over  W  code  chips 
for  each  phase  comparison,  then  this  synchronization  process  will  require  inte¬ 
gration  over  as  many  as  WN  chips. 

The  suggested  technique  for  the  synchronization  of  multi-state  maximal 
linear  codes  makes  use  of  their  relatively  high  and  regularly  positioned  auto 
correlation  sidelobes.  The  synchronization  procedure  uses  the  conventional 
search  technique  to  locate  one  of  the  set  of  equally  spaced  correlation  side¬ 
lobes.  Since  the  auto  correlation  sidelobes  are  spaced  M  chips  apart,  at 
most  M  relative  phase  shifts  of  the  two  (2)  codes  must  be  compared  in  order 


* 
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to  find  one  of  the  auto  correlation  sidelobes.  The  location  of  one  of  the 
auto  correlation  peaks  determines  the  phase  position  of  the  other  regularly 
spaced  correlation  peaks.  Each  of  these  phase  positions  is  then  tested  to 
determine  phase  synchronism  between  the  received  and  locally  generated  codes. 

Since  the  correlation  sidelobes  which  are  being  searched  out  in  this  pro¬ 
cedure  have  a  normalized  value  8  of  less  than  one,  the  integration  process 
must  be  taken  over  a  greater  number  of  code  chips  in  order  to  compensate  for 
the  fact  that  the  two  (2)  code  phases  being  compared  are  not  perfectly  corre¬ 
lated.  Thus  if  the  integration  is  performed  for  time  T  in  the  case  of  perfect 

correlation  (0  =  1),  then  the  integration  time  must  be  T,  =  in  order  to 

8 

obtain  an  equivalent  amount  of  energy  out  of  the  correlator.  The  following 
argument  established  this  fact. 

Suppose  that  in  the  case  of  perfect  correlation  the  voltage  out  of  the 

2 

correlator  is  V  or  equivalently  the  normalized  power  is  V  .  If  we  integrate 

p 

over  a  time  T,  we  obtain  an  energy  level  of  E  =  V  T.  If  the  normalized  corre¬ 
lation  function  0  is  less  than  one,  then  the  voltage  at  the  output  of  the 
correlator  is  0V  or  equivalently  the  normalized  power  is  (0V)  .  In  order  to 

obtain  the  same  signal  energy  at  the  output  of  the  correlator  as  in  the  pre- 

2  2 

vious  case,  we  must  integrate  over  a  time  T,  such  that  (0V)  Tj  =  V  T  or 


We  now  compare  the  acquisition  times  required  for  binary  and  multi-state 
maximal  linear  sequences.  We  assume  that  the  binary  sequence  is  acquired  by 
integrating  over  W  chips  in  each  relative  phase  position  so  that  the  acquisi¬ 
tion  process  will  require  integration  over  as  many  as  WN  chips,  where  W  is 
the  code  period.  To  acquire  a  multi-state  maximal  linear  sequence,  we  first 
find  one  of  the  correlation  peaks.  The  maximum  number  of  chips  over  which 
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we  must  integrate  is  chips.  The  maximum  number  ot  chips  over  which  we 

/  N  \  W 

must  integrate  to  find  the  mean  correlation  peak  is  (-w-  - 

VM/  (1-9) 


1  • 


Thus  the 


total  maximum  number  of  chips  over  which  we  must  integrate  is 
WM_  (h\  W 

02  \M/  (l-Q)2 

of  chips  WN  over  which  we  must  integrate  to  acquire  a  binary  maximal  linear 


The  ratio  of  this  number  of  chips  to  the  maximum  number 


sequence  is 


/M  \  J_  +  J_  1 
W  e2  "  (7^7 


I  »  fp  k  *  n  « 

For  maximal  2k-state  sequences  N  =  2k  n  -  1  and  M  =  z — - 


2  -  1 


formula  thus  becomes 


The  above 


a 


2k  n-  ijl  (i-e)2 


For  8-state  maximal  linear  sequences  we  have 


a 


Using  the  average  value  of  9  over  the  six  (6)  auto  correlation  sidelobes,  we 
have  0  =  -83 


a 


•21  + 


243 


,3  -n 


-  1 


This  function  is  plotted  in  Figure  7. 
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Figure  7 
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Comparison  of  Binary  and  Multistate  Synchronization  Times 
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3.8  Power  Soectral  Density  of  a  Maxima]  Linear  Sequence  over  GF(2*) 


In  this  section  we  determine  the  power  spectral  density  of  a  maximal 

k 

linear  sequence  over  GF(2  )  in  terms  of  the  complex  valued  Q-function  which 
has  been  defined  as  27ri 

Lk  ,  A  "j  .  . .  2k 


{  k  A  -j  . 

Qj<(j)  =  (2  -  1  -  j +jU)  +1  where  w  =  e 

The  power  spectral  density  S  of  a  maximal  linear  sequence  over  GF( 2k )  will  be 


S(o)  =  1 

F(Qk)(j ) 

S(j)  =  (P  +  1) 

2k 

0k  0  2  u  i  j  s 

t  v*>«  2l-' 

s=0 

where  F(Qk)(j)  = 

The  form  of  the  autocorrelation  function  and  power  spectral  density  of  a  maxi 
mal  linear  sequence  is  illustrated  in  Figure  8. 

IS  r~  /  r\  \  t  '  \ 


F(Qk)(j) 


F(Q. )  has  period  2K -  1 


2(2-1) 


3(2-1 ) 


Figure  8.  Power  Spectral  Density  of  Maximal  Linear  Sequence  over  GF(2K) 
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The  power  spectral  density  S  of  a  sequence  a  of  period  P  is  found  via 
the  Weiner  Khinchin  Theorem  by  taking  the  discrete  Fourier  transform  F  of 


its  autocorrelation  function  9  where  the  Fourier  transform  is  given  by 


■  j. 

S(a)(j)  =  F(9(a))(j)  =  ]P  e(a)(k)  e 


2tt  i  j  k 


In  this  section  we  show  that  if  a  is  a  maximal  linear  sequence  of  period 
P  =  2n'k-  1,  then 


S(a)(o)  =  1 


S(a)^j(2k  -  1)^  =  2kn  =  P  +1 

k(n-l)  . 

S(a) ( j )  =  2  F(Qk){j)  for  j  is  not  a  multiple  of  2-1 

We  first  make  the  following  computation  for  Qk 


Result: 


Let  Qk(j)  =  (2  -  1  -  j)  u>  +  j  w 


Then  Qk(j) 
j=0 


•(j+1) 


u  =  e 


Proof: 


T,  *J  5  H1- 

j=0 


k  k  k 

2-1  2  2-1 


X  -  X  -  X  +  1 
(1  -  X)2 


Differentiating  and  multiplying  by  X  we  have  the  identity 


k  2k  k  9k  i 


V*  jXJ  _  (2  -2m  -  (2  -1)X^  +  X 

U  U  -  x)? 
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Using  these  identities  with  X  =  w  we  obtain 


2k-2 

(i)  (2k - 1)  y 

j=0 


2(2k-l)  -  -  (2k-l)ui 

(1  -  u))2 


2k-2 

(2)  -  y  j“j 

j=0 


(2k-2)  -t-  (2k-l)cii~1  -  u> 

(1  -  o>)2 


(3) 


£ 


j=0 


(2k-2)u)  -  (2k-l)  +  J 

(1  -  u,)2 


(1)  +  (2)  +  (3)  =  1,  and  hence  the  result  follows. 

We  now  compute  S ( a ) 

P-1  _  2ri js 

Proof:  S(a) ( j )  =  F(0)(j)  =  ^  9(s)  e  P  P  =  2nk  -  1 

s=0 


Let  9'  be  the  autocorrelation  function  of  a  binary  maximal  linear  sequence 
of  period  2n  k-l.  It  is  well  known  that 


9'  (0) 
9*(s) 


2n"k  -  1  F ( 9 1 ) ( 0 ) 

and 

-1  for  s  t  0  F{ 9 ' )(j) 


1 

P  +  1  =  2kn  for  j  i 
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Since  0(s)  =  -1  for  s,  not  a  multiple  of  M  = 


2n-k  _  j 


2k  -  1 


L ,  we  have 


P  -1  2tr  f  j  S 

p 

s=o 


p-i 

2tt  i  j  s 

2-2 

23  e’^s^ e 

p  + 

23  D  +  0(sM)]e 

s=0 

S=1 

P-1 

2iri  js 

2k-2 

23  8’^s^ e 

p 

+  23 [i  +  e(sM)] 

s=0  s=l 


2tt  i  j  Ms 
P 


2th  s 
2k-  1 


Since  0(sM)  =  2k(n_1^  Qk(j)  -  1,  we  have 


•s. 


P-1 

23  0* (s) e 
s=0 


2»ijs  2k-2 

P  +  2k(n-1}  ^  Qk(s)e 

s=l 


2tt  j  j  s 
2k-l 


Since  Qk(0)  =  2k, 


we  have 


P-1 

23  O'OO  e 

s=0 


2tt  i  js 
P 


+  2 


k(n-l) 


2k-2 


Z  qk(s) 

s=0 


_  2nijs 
2k-l 


=  F(0')(j)  +  2k(rM)  F(Qfc)(j)  -  2k'  n 


i 

-  -< 

<  ; 
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Since  F(9‘ ) (0)  =  1  and  F<Qk)(Q)  =  2k,  we  have 

S(a) (0)  =  1 

F(6')(j)  =  P  +  1  =  2k  '  n  for  j  f  0.  Hence 

S(a)(j)  =  2k(n_1)  F(Qk)(j)  for  j  t  0 

2k-2 

•  Y  vs>  ■  2k 

s=0 

Thus  S(a)^j(2k-1)^  =  2kn  =  P  +  1. 

The  power  spectral  density  of  a  maximal  linear  sequence  of  period  63 
is  shown  in  Figure  9. 
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10  Log  [Power  Spectral]  of  Maximal  Sequence  over  GF(8)  Period  63 


APPENDIX  A 


Theorem: 


Proof : 


f (O)h(i) 


Let  4>  be  any  linear  mapping  of  the  field  GF(2°)  into  the  field 

k  k 

GF(2  ).  Let  f  be  any  polynomial  with  coefficients  in  GF(2  ),  i.e., 

f  c  GF(2^)[x].  Let  a  be  a  root  of  f  in  GF(2n).  The  sequence 

h(i)  =  if  (a  1)  is  a  member  of  V(f),  i.e.,  the  sequence  h  can  be 

generated  by  the  shif*  register  corresponding  to  the  polynomial  f. 


We  assert  that  the  sequence 

h(i)  =  f(a_i) 

can  be  generated  by  the  shift  register  corresponding  to  the  polyno¬ 
mial  f  (x)  where  f(a)  =  0.  The  sequences  h  generated  by  the  poly¬ 
nomial  f(x)  are  those  which  satisfy  the  recursion  relation: 

=  f  (l)h(i-l)  +  f  (2)h(i'-2)  +  .  .  .  f  (n)h(i-n)  for  all  i  >  n 

Thus,  to  show  that  the  sequence  h  can  be  generated  by  the  polynomial 
f,  we  must  show  that  is  satisfies  the  recursion  relation. 


f(a)  =  0 


hence  f(0)  +  f(l)a  +  f ( 2 ) a.2  +  .  .  .  f(n)c.n  =  0 


-1 


multiplying  by  a  we  have 


f (0)a_i  +  f(l)a"(l_1)  +  f(2)cT(1"2)  +  .  .  .  f (n) a" (l“n)  =  0 


♦  (0) 


^f (0)a-i  +  f(l)a~(1_1)  +  .  .  .  f (n)a~(i~n)j  = 

(a-1)  +  f(l)f  +  .  .  .  f(n)^a'(1  "n)^  -  0 


f(0H 

f(0)h(i)  +  f(l)h(i-l)  +  .  .  .  f (n)h(i-n) 


A-l 


Theorem 


k  k 

(Converse)  Given  any  2  -state  sequence  h  with  values  in  CF(2  ) 

which  has  been  generated  by  a  shift  register  corresponding  to  the 

irreducible  polynomial  f(x)  of  degree  n  and  having  coefficients  in 

R  n  k 

GF(2  ),  then  there  is  a  linear  mapping  $  of  GF(2  )  into  GF(2  )  such 

that  h(i)  *  $(a  1)  where  a  is  a  root  of  f(x). 


A-2 


APPENDIX  B 


In  this  appendix  we  review  the  properties  of  the  projection  mapping  [l] . 

Result:  For  all  Y  e  GF(22)  we  have  P2(y)  -  Pj^B  *Y^ 

P3Cr)  =  P^BY) 

Proof:  Y  =  a(0)  +  a(l)B  +  a(2)B2 

BY  =  a (0)B  +  a(l)B  +  a(2)(B2  +  1) 

=  a  (2)  +  a  (0) B  +  (a ( 1 )  +  a(2)j  B2 
6_1Y  =  a(0)B6  +  a ( 1 )  +  a(2)B 

=  a (0) (B  +  B2)  +  a ( 1 )  +  a(2)B 

=  a ( 1 )  +  a (0)  +  a (2)  B  +  a(0)B2 

Pj (B_1 Y)  =  a( 1 )  =  P2(Y) 

Result:  f(x)  e  GF(23)[x]  we  have  P^  (f  (x))  =  Pj(p  *f(x)) 

P3(f(x))  =  P  j  (b  f  (x)) 

Proof:  Reference  Gold  Notes  (Page  80) 

P^fU)) 

Pj(b  I  f  (i)x1) 

Pj(l(6f  (i)V) 


EP^efU))  x1 

definition  of  P 

lP3(f(i))  x1 

previous  result 

P3(lf(i)x^ 

P3(f(x)) 

[l]  Study  of  Multi-State  PN  Sequences  and  Their  Application  to  Communica¬ 
tion  Systems.  Contract  N00014-75-C-1040  -  23  July  1976. 
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Result : 


Let  a  = 


Proof  : 


Result 


Proof  : 


g(x) 
f  (x) 


be  a  rationalized  representation  of  a  multi-state 


sequence  where  g(x)  e  GF(2  ) [x]  and  f(x)  e  GF(2)[x].  The  represen- 

p  (  x  ) 

tation  of  a  =  ■  is  the  interleaving  of  k  binary  sequences  is 


given  by: 


a  = 


P1(g(x))  P2(g(x)) 


f  (x) 


f  (x) 


Xs(x)) 

f  (x) 


Reference  Gold  Report  (Page  83) 


a 

g 

f 

a  •  f 

=  g 

P  (a  •  f) 

=  Pi(g) 

f  •  P.(a) 

«  p.(g) 

P.(a) 

pi(g) 

f 

r 

1 

P. 

1 


„k  g(x) 
J  '  f(x) 


f  e  GF( 2 ) 

see  result  in  Gold  Notes  (Page  81) 


=  P. 


Bkg(x) 


f(x) 


ek  •  a 


(B  )  (a) f (x) 
(Bk)(a)f(x)J 
f(x)Pi^6k)  (aj) 


Pil  U 


kUg(x)> 


g(x) 
f(x) 

„k\  I SOO 

f  3 


B 


( x  V 


(B  )  g(x) 

Pi 

Pi(!ek)  g(x)) 

Pj  ((fk )  g(x)) 
f  (x) 
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APPENDIX  C 


AUTOCORRELATION  OF  MAXIMAL  LINEAR  SEQUENCES  OF  PERIOD  63  OVER  GF(8) 
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4  2  5  31 

Sequence  Polynomial  ax  +  a  x  +  1,  Root  a 


-  RtAL  -  . IH  AG  INAPT"  — 

I  INI.  (SAP  I  2)/2  INI,  (SPTvf  2)/2 
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17  -1 
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20  -1 
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22  -1 

23  -1 

20  -1 

2P  -1 
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30  -1 
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6  2  6  22 
Sequence  Polynomial  ax  +  a  x  +  1,  Root  a 


— 

REAL  - 

-  IMAGINARY 

— 

I 

IWT. 

(SORT  2 !/7 

INT.  (SQR! 

2)/2 
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4  ?  4 

Sequence  Polynomial  a  x  +  ax  +  1 


- REAL  -  - IMAPINARY - 

I  INI.  (SORT  2)/2  INI.  (SORT  2)/2 
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